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Abstract: An approximate homotopy symmetry method for nonhnear problems is pro- 
posed and apphed to the six-order boussinesq equation which arises from fluid dynamics. 
We summarize the general formulas for similarity reduction solutions and similarity reduc- 
tion equations of different orders, educing the related homotopy series solutions. Zero-order 
similarity reduction equations are equivalent to Painleve IV type equation or Weierstrass 
elliptic equation. Higher order similarity solutions can be obtained by solving linear vari- 
able coefficients ordinary differential equations. The convergence region of homotopy series 
solutions can be adjusted by the auxiliary parameter. Series solutions and similarity re- 
duction equations from approximate symmetry method can be retrieved from approximate 
homotopy symmetry method. 
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1 Introduction 

Nonlinear phenomena arise in many aspects of science and engineering. Lie group theory [1-4] 
provides remarkable techniques in effectively studying nonlinear problems such as exploring sim- 
ilarity reduction of partial differential equations. It should be noted that approximate solutions 
also contribute to understanding the essence of nonlinearity. Perturbation theory [5-7] was 
consequently developed and it plays an essential role in nonlinear science, especially in finding 
approximate analytical solutions to perturbed partial differential equations. 
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Combined with Lie group theory, perturbation theory gives rise to two distinct approximate 
symmetry methods. For the first method due to Baikov, et al [8,9], symmetry group generators 
are generalized to perturbation forms. The second method proposed by Fushchich, et al [10] 
is based on the perturbation series for the dependent variables which decomposes the original 
equation into a system of equations. Approximate symmetry of the original equation boils down 
to exact symmetry of such system of equations resulted from perturbation. The comparisons in 
Refs. [11, 12] show superiority of the second method to the first one. 

Aside from perturbation theory, some nonperturbative techniques, such as the artificial small 
parameter method [13], the (5-expansion method [14] and the Adomian's decomposition method 
[15] are also of significance when perturbation quantities are not involved in many problems. 

Liao [16] developed a kind of analytic technique, the homotopy analysis method, and solved 
some problems successfully [17-19]. Recently, the homotopy analysis method was further im- 
proved in Refs. [20-22]. The homotopy analysis method is based on homotopy conception in 
topology. This method is suitable for problems that contain no small parameters. Further- 
more, the series solutions obtained from the perturbation method, the artificial small parameter 
method, the 5-expansion method and the Adomian's decomposition method can also be retrieved 
by the homotopy analysis method. 

Perturbation techniques are only confined to weak perturbation problems. For strongly per- 
turbed nonlinear systems, we propose approximate homotopy symmetry method to study possi- 
ble analytic series solutions. The six-order boussinesq equation is used as an example to illustrate 
the effectiveness of the homotopy symmetry method. 

2 Basic notions 

For a nonlinear partial differential equation 



where A is a nonlinear operator, u = u{x, t) is an undetermined function, and {x, t} are 
independent variables, we introduce a homotopy model 



A{u) = A{X, t, Ux, Ut, Uxx, Uxt, ■■■) =0, 



(1) 



niu,q)=0, 



(2) 
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with q G [0, 1] an embedding homotopy parameter. The above homotopy model has the property 



n{u,i) =A{u), n{u,o) =no{u), 

where TCoiu) is a differential equation of which the solutions can be easily obtained. 
We make an ansatz that the homotopy model ([2]) has the homotopy series solution 



(3) 



where Ui solves the system 



E 

i=0 



Tio{uo) = 0, 

n'o{uo)ui + Fi{uo) = 0, 

'H'q{uq)u2 + -^2(^0, Ul) = 0, 



TiQ{uo)ui + Fi{uo, Ul, 



Ui-i) = 0, 



(4) 



(5) 
(6) 
(7) 

(8) 



in which the operator Ti.Q{uo) is defined as 



d 



e=0 



for arbitrary function f{x,t), and all Fj = Fj(no, ui, • • • , Uj-i) satisfy 



1 & 



Ukq 



(i = l, 2, •••). 



q=0 



Then, the solutions of the original nonlinear system ([T]) read 

00 



1=0 



(9) 



(10) 



Now, we are concerned about constructing approximate group invariant solutions of the 
homotopy model ([2]). First, we introduce the following definitions: 

Definition 1. Symmetry (or exact symmetry). A symmetry, cr, of the nonlinear equation ([T]) 
is defined as a solution of its linearized equation 



d 

A'{u)a = —A[u + ea] 



0, 



(11) 



£=0 



which means that Eq. ([T]) is form invariant under the transformation u ^ u + ea with an 
infinitesimal parameter e. 

Definition 2. Homotopy symmetry. A homotopy symmetry, ag, of the nonUnear system ([1]) is 
an exact symmetry of the related homotopy model i.e., a solution to the linearized equation 
of @ 

n'{u,q)aq = 0. (12) 

Definition 3. Approximate homotopy symmetry. The A;th order approximate homotopy sym- 
metry, (Tg^fc = ((To, (Ti, C72, • • • , crfc) ) of the nonlinear system ([1]) is an exact symmetry of the 
system of the first (k + 1) approximate equations, i.e., the solution of the following linearized 
system 

K{uo)ao = 0, (13) 

fc-i 

n'o{uo)ak + iKiuo))' aouk + (Fk)'^^ = 0, (j = 1, 2, k), (14) 

j=0 

where the operator {Ff^)'^, is defined as 

d 

i^k)L = -K-Fk{uo, Ml, Uj+eaj, uj+i, n^^i) 

e=0 

For simplicity later, the following simple homotopy model is exclusively taken 

{1 - q)no{u) + qXA{u) = 0, (15) 

with A 7^ an auxiliary parameter. It is easily seen that Eq. (|15p varies asymptotically from 
Tio^u) = to Eq. ([I]) as (7 goes gradually from to 1. When Tio is fixed as a linear operator, 
the homotopy model ([2]) is just the usual one applied in Refs. [16-22]. 

From the above process, we see that the approximate homotopy symmetry method is an 
integration of the homotopy concept, perturbation analysis and symmetry method. 

3 Approximate homotopy symmetry method to the six-order 
boussinesq equation 

The illposed Boussinesq equation 

Utt = {u + U^)xx + Uxxxx, (16) 



describes propagation of long waves in shallow water under gravity [23], in one-dimensional 
nonlinear lattices and in nonlinear strings [24]. Filtering and regularization techniques was 
applied to Eq. (fT6]) in Ref. [25] to introduce the singularly perturbed (sixth-order) Boussinesq 
equation 

Vtt = Vxx + + 

^xxxx ~l~ £ ^XXXXXXl iX^^ 

where e is a small parameter. In Ref. [26], double-series perturbation analysis was utilized to 
recover Eq. (fT7|) . 

For formal brevity, we rewrite Eq. (jl7|) as 

^tt ~\~ '^xx ~l~ i,U )xx ~l~ Uxxxx — IJ^'Uxxxxxxi (-'-^) 

where u is a function with respect to x and t, ^ is an arbitrary parameter. Assuming 'Ho{u) to 
be the left hand side of Eq. (jlSp and replacing A by 1 — for concision of the results, we change 
Eq. into 

(1 - q6){uu + Uxx + {u'^)xx + Uxxxx) - qfJ-C^ - 0)Uxxxxxx = 0. (19) 

It is easily seen that Eq. (fT9]l is just the boussinesq equation [27,28] when g = 0. 

Substituting Eq. into the above equation and matching the coefficients of different powers 
of q yield the following system of partial differential equations (/c-order approximate equation) 

k k-l 

Uk,tt~\~Uk^xx~\~'^ ^ '] {U'k~i,xU'i^x~^Uif—iUi^xx)~^Ui^ xxxx~^t^i^ 1) ^ ^ ^ Ui^xxxxxx — (^ — 0, 1, " " " ) 

i=0 i=0 

(20) 

with M_i = 0. 

We investigate similarity reduction of Eq. (j20p by Lie symmetry method [29] . The linearized 
equations for Eq. ([20]) are 

k 

^k,tt ~l~ ^k,xx ~l~ 2 ^ '] i^k—i,xU,i^x ~l~ ^fc— ~l~ ^k—i^i^xx ~\~ Up^^iGi xx) 
i=0 

k-l 

+(^k,xxxx + /^(^ ~ 1) ^ ^(yi,xxxxxx = ^, (^ = 0, 1, •••) (21) 

1=0 

where cr^ are functions of x and t with cj-_i =0. Eq. (j2ip means that Eq. (j20p is invariant under 
the transformation Uk ^ Uk + ecTk, (k = 0, 1, • • • ) with e an infinitesimal parameter. 



5 



The symmetry transformations 

ak = Xuk,^ + Tuk,t-Uk, {k = 0, 1, •••), (22) 

where X, T and are functions with respect to x, t and u/^, (k = 0, 1, • • • ), conform to Eq. 
([2T]) under Eq. ([20]) . We consider finite equations in Eqs. ([20]) . (pT|) and ([22|) to summarize 
general formulas for similarity reduction solutions and similarity reduction equations. 

Confining the maximum of to 2 in Eqs. (|2U|) . (|21|) and (|22|) . we see that the independent 
variables of X, T, Uq, Ui and U2 are accordingly restricted to x, t, uq, ui and U2. More than 
2000 determining equations are obtained by inserting Eq. ()22p into Eq. (j2ip . eliminating Uo.tt; 
ui^tf and ?i2,jt in terms of Eq. (|20|) and vanishing coefficients of different partial derivatives of 
liQ, ui and ti2. 

To solve the determining equations, we first extract the simplest equations for T 

from which we have T = T{t). Considering this condition, we select the simplest equations for 
X 

Xt = XuQ = Xui = Xu2 = 0, 
with the solution X = X{x). In this case, we get the simplest equations for C/qi Ui and U2 

f^2,tMo — f^2,i?ii — ^2, Mono — ^2,tiowi — f-^2,uoM2 — ^2,miui — U2^u\U2 — U2^U2U2 — 0) 

which imply 

UQ = Fi{x,t)uo + F2{x,t), 
Ui = F3{x)uo + F^ix, t)ui + F5{x)u2 + Fe{x, t), 
U2 = Ft{x)uo + F8{x)ui + FQ{x,t)u2 + Fio{x,t), 
where the undetermined functions satisfy 

F3 = F^ = Fq = Ff = Fio = Fi^x = Fi^t = F2^x = -^2,** = = -Rl,t = -^8,a; = Fg^x = Fg^t 
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— — Ttt — 0, Tt — 2Xx, 2Xx — 2Tt + Fi, — Tt + F2, 2Xx — 2Tt + 2F4 — Fq, 

6Xx = 2Tt + F4 - Fg, 6Xx = 2Tt + Fi - F4, 69X^ = 2eTt + e{Fi - Fg) - Fg. 

The solutions to the determining equations are finahy obtained by solving the above system as 
follows 

X = ^Cix + C3, T = Cit + C2, Uo = -CiUo-^Ci, Ui = -2CiUi, U2 = Ci{eui-3u2), (23) 

where Ci, C2 and C3 are arbitrary constants. 

In the same way, limiting the maximum of A; to 3 in Eqs. (I20p . (I2ip and ()22p . we execute 
similar computation and obtain 

X = ^Cix + C3, T = Cit + C2, f/o = -Cmo - \ci, 

Ui = -2CiUi, U2 = CiiOui - 3u2), Us = 2Ci{eu2 - 2^x3), (24) 

where Ci, C2 and C3 are arbitrary constants. 

Enlarge the domain of k by degrees and repeat similar procedures , we discover the formal 
coherence of X, T and Uk {k = 0, 1, • • • ) , i.e., 

X = ^Cix+Cs, T = Cit+C2, Uk = Ci[{k-l)9uk-i-{k+l)uk]-^Ci5k,o, (A: = 0, 1, • • • ) (25) 

where C^, {k = 0, 1, •••) are arbitrary constants. The notation 5kfi satisfying 5o,o = 1 and 
^k,o = (A; 7^ 0) is adopted in the following text. Similarity solutions to Eq. (120]) can be obtained 
by solving the characteristic equations 

dx dt duQ dt dui dt duk dt 

Y ~ ~ Ih ~ "' ' ~ "' ^ ' 

which are distinguished in two subcases. 

3.1 Homotopy symmetry reduction of Painleve IV type 

When Ci 7^ 0, without loss of generality, we rewrite the constants C2 and C3 as C1C2 and 
^CiCs and change Eq. (f25]) to 



X = ^Ci{x + C3), T = Ci{t + C2), Uk = Ci[{k-l)euk-i-{k + l)uk]-^Ci5k,o, {k = 0, 1, 



(27) 
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From the first three equations in Eq. ()26p . we get the invariants 

''^•" = « = 7R%' '''' 

Io{x,t,uo) = Po = uot + C2U0 + ^t, (29) 

Ii{x,t,ui) = Pi = {t + C2fui. (30) 

Viewing Pq and Pi as functions of ^, we have 

^0 = ITTTTTTT' (31) 



2(t + C2) 

(i + C2)2- 



Similarly, we get other similarity solutions 



"3 = , ^ X. + , ^ X. + , (34) 



(t + c72)4 (t + c2)4 (t + a 



2f' 



which conform to the general expression 

2^0(0- ^ ^^.fc-n(^im^ a n 1 ^ C^^^ 

2(t + g2) + ^ ^ ^ ^T^Tcy^' ^^ = '' 

with the similarity variable ^ = 

From Eq. (jH), we get the series solution to Eq. ([19]) 

PkiOQ' 



u 



2(* + C2) (t + C2)(t + C2-etg) 



and when further setting g = 1, we have 



2(* + C2) (t + C2)(t + C2-0t)'=' 

which is a homotopy series solution to the six-order boussinesq equation. 

The determination of similarity reduction equations depends on finite equations in Eqs. (j20p 
and psp . It should be emphasized that all the previous similarity reduction equations should 
be considered to remove the terms Po,^^^^) -f'l,????' ' " " ' rather than Pg^^^^^^^, Pi,^^^^^^, 
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• • • , Pfc_i^gggggg when we eliminate Uk in the kth. order approximate equation (j20p in terms of the 
similarity solutions (|35|) . We sum up the general formula for the similarity reduction equations 

-{k - i)C2eiPk-u + + i){k + 2)Pk - 2{k + i){k - i)C2ePk-i + {k- i){k - 2){C2efPk-2 

k-l 

+C2Sk,o + ^K0-l)iC2e)''-'-'P^,^^^^^^ = 0, (fc = 0, 1, •••) (38) 
with P„2 = P^i = 0. 

When k = 0, Eq. (j38p is equivalent to the Painleve IV type equation. When k > 0, specific 
forms of Eq. (|38p depend on the solutions Pq, Pi, ■ ■ ■ , Pk-i and we can rearrange the terms in 
Eq. JMD as 

^k,£,iii + '^{PkPo,(^ + PoPk,^() + 4Pfc,5Po,C + + C2)Pk,S,i 

+^(4fc + 7)CPfc,5 + (A; + + 2)Pfc = (A: = 1, 2, • • • ) (38') 

where fk{Cj is a function of {Pq; Pi-, " ■ ■ > -Pfc-i} 

fkiO = {k- l)C2eiPk-U + 2{k + l){k - l)C2ePk^i -2Y,{Pk-^P^,ii + Pk-i,^P^,^ 

1=1 

k-l 

- 1)(A: - 2)(C20)2Pfc_2 -^^0- l)iC2e)'-'-'^,iaai. 

1=0 

Eq. (38') is actually a fourth order linear variable coefficients ordinary differential equation of 
Pfc when Pq, Pi, • • • , Pk~i are known. 

Remark: Taking 6 = and making the transformations Pk{0 — f^^QkiO — ' ' ' ); 

change the similarity reduction equations (p8]) and the homotopy series solution ([37|) to 



£ 1 

Qk,!;m + 2 X](*5fc-i<5i,cc + Qk-i,^Qi,^) + ( J + C'2)Qfe,e5 + ^ (4A; + 7)eQfc,^ 

+{k + l)(/c + 2)Qfc + C2Sk,o - Qk-imm = 0' = 0' 1> • • • ) (39) 



and 

2(t + C2)+^Jt + C2)'=+i' ^^^^ 
respectively. It is easily seen that Eqs. (fTSj) . ([39]) and (flO]) coincide with Eqs. (4), (21) and (20) 
in [30] respectively. 



3.2 Homotopy symmetry reduction of the traveling wave form 

When Ci = 0, from the characteristic equations (j26p . we get the similarity solutions 

uk = pm, (fc = o, 1, •••) (41) 

where ^ = t — ^x. We take an equivalent travelling wave form = x + ct for the similarity 
variable with c an arbitrary velocity constant. 

From Eqs. (jU and (fTO|) . we obtain a series solution to the homotopy model ([19]) 

oo 

u = j2i''m), (42) 

k=0 

and a homotopy series solution to the six-order boussinesq equation 

oo 



oo 

U = 

k=0 



By substituting the similarity solutions (j^Tj) to approximate equations ([20ll , we get the simi- 
larity reduction equations 

k k-1 

PkmiHc^+'^)Pk,ii+2^{Pk-iPi,i^+Pk-i,iPi,i)+K0^ (fc = o, 1, •••) 

1=0 j=0 

(44) 

which are equivalent to 

Pk,ii + (c' + l)Pk + 2(2 - 5k,o)PkPo = 9ki0, ik = 0, 1, • • • ) (44') 

where 

k-l k-l 

gkiO = -2 ^ - fi{9 -1)Y^ e^-'-'Pi,^^^^ - aki - bk, 

i=l i=0 

with ttk and 6^ arbitrary integral constants. Eq. (44') is a second order linear variable coefficients 
ordinary differential equation of Pk- 

Remark: Taking 6 = and making the transformation Pk{^) = fJ'^QkiO {k = 0, 1, •••), we 
reduce the similarity reduction equations (01]) and the homotopy series solution to 

k 

Qkm^ + i^'^ + '^)Qk,iiii + '2Y.iQk-iQi,^^ + Qk-i,^Qi,^)-Qk-imi;i;i; = ^^ (A: = 0, 1, •••) (45) 

1=0 

and 



oo 

U = 

k=0 
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respectively, with Q-i = 0. It is easily seen that Eqs. (fT8|) and (|i6|) are identical to Eqs. (4) 
and (24) in [30] respectively. Integrating Eq. (|45|) twice, we obtain Eq. (25) in [30]. 

When k = 0, Eq. (j44p is equivalent to the Weierstrass elliptic equation and has the following 
solutions 

Po(0 = - i(c2 + 1) - tanh2(a^ + b), (47) 

Po(0 = 2a2(l - 2m?) - ^(1 + c^) + 6a2m2cn2(ae + b), (48) 

^'o(e) = -^(l + c^) -6a2p« + 6,m2,m3), (49) 

where a, 6, mi, m2 and ms are arbitrary constants. 

When A: > and oq = 0, Eq. (44') can be integrated out one after another 



ek + j Po,| (cfe + J Po,i9kd^^ de , 



(50) 



where Cfc and are arbitrary integral constants. 

In the rest of this section, we discuss the solutions of similarity reduction equations (j44[) 
from hyperbolic tangent function solution (j47p . We suppose that Eq. (j44p have the hyperbolic 
tangent function solutions 

^fc(e) = J^4,itanh*C, ik = 0, 1, •••) (51) 

1=0 

where all d^^i are constants to be determined. By balancing the highest powers of tanh^ from 
Pfc_i^^^^^^^ and -Pfc,55 in Eq. dM]), we have n^.i + 6 = + 4, with no = 2 in Eq. ([^7]) . leading 
to Uk = 2{k + 1). 

The kth similarity reduction equation in Eq. ()44p contains Pq, Pi, • • • , P^) so that Pq, Pi, 
■ ■ ■ , Pk can be solved one after another starting from Eq. (j47p . For the kih. equation in Eq. 
(j44p ■ Pq) -Pi) ■ ■ ■ ) Pk-i fire known and inserted into Eq. (j44p together with Eq. (j5ip . Matching 
the coefficients of different powers of tanh^, we get a system of algebraic equations with respect 
to dkfi, • • • , dfe,2{A;+i) from which we can construct Eq. ([5T]) . 

We list dk^i in Eq. ([5T]) up to A; = 3 
A; = 1 : 



di^o = 22^a^{0 — 1), di^i = pi, (ii^2 = — 12O/^a'^(0 — 1), di^a = — pi, (ii^4 = 90/ua' 



4/ 
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2 

d2,o = l8Ofi'^a%0 - if + 22iia^e{9 - 1) - ^,d2,3 = -SOna^piiO - 1) - P2, 
d2,i =P2, d2,2 = -2970/^2 a^(0 - 1)2 _ l2O^a^e(0 - 1) + ^2,5 = 3O^aV(0 - 1), 

2 

^2,4 = 558O//2a6(0 - 1)2 + gOfja'^e^e - 1) - ^^2,6 = -279O/i2a*^(0 - 1)2, 



A; = 3 : 



,3^8//) 1^3 , Q«n,,2„6/)/'/) i \2 , oo,,^4/)2m 1^ i ^..„2/fl i\ ^'lP2 



dso = 6030^^ a**(6l - 1)^ + 360//^a''6'(6l - 1)^ + 22/ua*6l^(6l - 1) + -^K,. 

6 12a^ 

d3 2 = -llSSOO/i^a^fe - if - 5M0n'^a^9(e - if - 120^ia^e\9 - 1) + —fipjie - 1) + 

12 3a^ 

3 

d3,3 = -1260/^2 a - 1)' - 30f,a\p^e + P2){e - 1) -^3 + 
^3,4 = 361080^^0^(0 - if + 11160fi^a^9{e - if + gO/ua^^^^^ _ _ i^^^a^^ _ P1P2 



2 ^^^^ ^ 4a2 



3 

(ig 5 = 2655/1^0^(6' - 1)^ + 30/io^(pi6l + p2)(6' - 1) 



72a4^ 



25 

d36 = -397296/i^ 0^(6* - 1)^ - 5580fj,^a^e{e - if + —npji9 - 1), 
d^^^ = p3, ^3,7 = -1395 fi'^a'^pi{9 - if, d3,8 = 148986//^ 0^(6* - l)^ 

where pi, p2 and ^3 are arbitrary constants. 

We specify the parameters by a = 1, 6 = 0, c = 1, = 0.1, pi = 0, p2 = 0, P3 = 0, p4 = 
and p5 = 0. To see the relationship between the convergence of the homotopy series solutions 
and the parameter 9, we choose four values 0.4, 0.8, 1.2 and 1.6 for 9 and display the plots of 
l-ffc(OI = 1) 2, 3) for —5 < X < 3 and t = 1 in Figure 1 where the dotted line, the dashed line 
and the solid line represent |Pi(.^)|, |-P2('?)I ^^id |-P3('?)| respectively. The necessary condition for 
the convergence of homotopy series solutions ()43p is lim^^oo |-Pa;(OI — and the convergence 
regions correspond to lower solid lines in Figure 1. 

For (a) and (b) in Figure 1, the convergence region of the homotopy series solutions (j43|) for 
9 = 0.8 is wider than the convergence region for 9 = 0.4. The homotopy series solutions (j43p 
corresponding to (c) and (d) in Figure 1 are divergent, with the homotopy series solutions for 
9 = 1.6 diverging faster than the homotopy series solutions for 9 = 1.2. 

Remark: When 9 < 1 {X > 0), the convergence regions of the homotopy series solutions (j43p 
grow wider provided that — > 1 (A — > 0). When 6* > 1 (A < 0), the homotopy series solutions 

12 



(j43p are divergent and the larger the values of 9 (the smaller the values of A), the faster the 
homotopy series solutions diverge. 

4 Summary and discussion 

In the framework of approximate homotopy symmetry method, we investigated the six-order 
boussinesq equation and summarized the similarity reduction solutions and similarity reduction 
equations for approximate equations of different orders. The homotopy series solutions to the 
six-order boussinesq equation were derived. 

Zero-order similarity reduction equations are equivalent to Painleve IV type equation and 
Weierstrass elliptic equation, fc-order similarity reduction equations are linear variable coef- 
ficients ordinary differential equations of Pki^,) which depend on particular solutions of the 
previous similarity reduction equations from zero-order to {k — l)-order. 

For homotopy symmetry reduction of the traveling wave form, we constructed hyperbolic 
tangent function solutions to /c-order similarity reduction equations {k = 1, 2, 3). |-Pfc('?)| were 
plotted for different values of A. 

The auxiliary parameter A dominates the convergence regions of the homotopy series solu- 
tions. For A > 0, the convergence regions get wider when A ^ 0. For A < 0, the homotopy series 
solutions are divergent and smaller values of A correspond to faster diverging homotopy series 
solutions. When A — > 0, Eq. (fT9l) tends to the boussinesq equation. Accordingly, the homotopy 
series solutions to the six-order boussinesq equation behave more and more like the solutions 
to the boussinesq equation. This may account for the role of A in controlling the convergence 
regions of homotopy series solutions. 

The approximate symmetry method is valid only if the nonlinear problems contain small 
parameters. In contrast, the auxiliary parameter A for the approximate homotopy symmetry 
method is artificial introduced. In other words, the approximate homotopy symmetry method 
is applicable to those nonlinear problems that contain no parameters at all. For the six-order 
boussinesq equation, the convergence of all approximate symmetry series solutions in Ref. [30] 
is influenced by the parameter fi. However, the parameter fj, is not explicitly included in the 
homotopy series solutions (j37p and (j43p . showing that the convergence of the homotopy series 
solutions is not directly affected by large parameters. 
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Moreover, the series solutions from approximate symmetry method can also be retrieved by 
approximate homotopy symmetry method. Therefore, the approximate homotopy symmetry 
method is superior to the approximate symmetry method and provides an effective analytical 
tool in constructing series solutions for many nonlinear problems. 
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